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Abstract 

We use a diquark model of the nucleon to calculate the electromagnetic form factors 
of the nucleon described as a scalar and axialvector diquark bound state. We provide 
an analysis of the zero-mode contribution in the diquark model. We find there are 
zero-mode contributions to the form factors arising from the instantaneous part of the 
quark propagator, which cannot be neglected compared with the valence contribution 
but can be removed by the choice of wave function. We also find that the charge 
and magnetic radii and magnetic moment of the proton can be reproduced, while the 
magnetic moment of the neutron is too small. The dipole shape of the form factors, 
g m(Q 2 )/r p and G^^Q 2 ) / /j, n , can be reproduced. The ratio [iG p E /G p M decreases with 
Q 2 , but too fast. 
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1 Introduction 


Methods using relativistic formalisms directly without Foldy-Wouthuysen reduction have 
been suggested HJ000 to include relativistic effects to all orders, while avoiding truncated 
power expansions in v/c or p/m. Among these relativistic formalisms, light-cone quark 
models are widely used. There are several works where the electroweak properties of mesons 
have been studied using light-cone quark models 0 0 0 El 13 DU- Applications of the 
Bakamjian-Thomas construction of the relativistic few-body system 0111 HD EE! to extend 
nonrelativistic wave functions to the relativistic domain are close to the naive quark model. 
Weber et al. nain give relativistic spin-flavor wave functions of a Lorentz covariant quark 
model. Araujo et. al. £0. connect these wave functions with the triangle Feyrnan diagram 
and give an effective Lagrangian of the nucleon, which corresponds to a field-theoretical 
formulation of the form factors. 

From the Eqs. (9,10) in Ref. [13f (KW, for short. See also Eq. (3) below.), we see that the 
nucleon wave function is split into two parts as (qq)q. At this point, the diquark model can 
be introduced naturally. The diquark model has been widely used to study the property of 
baryons. For example, nucleon properties are studied through the Bethe-Salpeter equation 
with scalar and axialvector diquark correlations M- Kroll et. al. study the electromagnetic 
form factors of the proton in the time-like region and two-photon annihilation into proton- 
antiproton in the infinite momentum frame mi. 

In all investigations of baryon properties 00111111211131111113, zero modes have not 
been considered. The zero-mode contributions can be interpreted as residues of virtual pair 
creation processes in the q + {= q° + q 3 ) —> 0 limit of the virtual photon momentum HM3. 
In the absence of zero-mode contributions, the hadron form factor can be obtained in a 
straightforward way by taking into account only the valence quark contributions. In an 
effort to clarify this issue, Jaus ra and Choi et. al. |2T independently investigated the 
electroweak form factors of spin-1 mesons in the past few years. Choi et. al. find that 
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the zero-mode contributions in principle depend on the form of the vector-meson vertex 
= 7 m — (P v — 2k)^/D] the form factor f(q 2 ) is found to be free from a zero-mode 
contribution if the denominator D contains a term proportional to the light-front energy 
( k~) n with a power n > 0. Analogs of T^ occur for the spinor parts of the nucleon current 
matrix element, as we shall discuss below. 

In this paper we calculate the electromagnetic form factors of the nucleon, taken as a 
scalar and axialvector diquark bound state, and investigate the zero-mode contributions in 
the quark-diquark model. 

In the following section, the nucleon wave function in the diquark model and an effective 
quark-diquark-nucleon effective Lagrangian are given. Then the extraction of the nucleon 
form factors from the current matrix elements are given in Section 3. The discussion of the 
zero mode is given in Section 4. The numerical results for the form factors in the frame 
where the virtual photon momentum q + = 0 are given in Section 5. The summary and 
discussion are given in the last section. 

2 Wave Function of the Nucleon 

In the static case the nucleon wave function with SU(6) symmetry can be written as 

-j={x P & P + X\® a)^o, (1) 

where Xp^p corresponds to the scalar-isoscalar component of the nucleon and to the 

vector-isovector quark pair part. This wave function can be rewritten in quark-diquark form 
as 

|pl) = —=[{2V^d 1 — \/2V^d} — VzVio + V^v}) sin a + 3Sqou^ cosa], (2) 
V18 

where the V are axial-vector diquark amplitudes and S is the scalar diquark amplitude. The 
explicit form of the V and S is given in Ref. [22] • If we set a = 7 t / 4 , then Eq. ([2]) reduces to 
Eq. (fTj) . that is, the permutation symmetry is restored. 
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Following KW jQIj, the two parts of the wave function in the nonstatic case become 

Xp®p -> N s (qn^Cqj) ( q 3 N x ) ; 

-> N v (qrfriTzCqZ') • {%XpTl ^ N \), (3) 

where a sum over quark permutations is implied. Here C = ? 7 2 7 ° is the charge conjugation 
matrix, Ns, Ny are the normalization constants of the diquark wave functions. Also N\ = 
[P]u\ with u x (P) the Dirac-spinor of the nucleon and [P] denotes 7 ■ P + tun with mjv the 
nucleon mass, g* = [P]ui with tq the quark spinors; a factor [P] with each quark spinor 
appears with the Melosh rotation of the nonrclativistic spin wave function to the light cone. 
This is derived in detail in Sect. IV.A of Ref. [25. . 

Here we recognize and gi 7 At rzr 2 CgJ as the scalar and axialvector diquarks 

from the NJL model [21] • Thus, the nucleon wave function can be written as (see the initial 
quark-diquark vertex in Fig. |TJ) 

\N) = fs<Ps{x, k ± )u(k)O s u x (P) + fv<Pv{x, k ± )e l a u(k)Oy^=ux(P), (4) 

= 1, 0“ = ( 7 ° + P“/M) 7 5 , (5) 

where ips , and tpy are the scalar and axialvector radial diquark wave functions and e a is the 
axialvector diquark polarization vector, and k M is the quark four-momentum. If coupling 
constants fs and fy is free completely, there will be a common factor, which will give a 
non-one charge of proton, so they must be constrained to give correct charge of proton. 
Our wave function is similar to that of Kroll et al. HZI and its basic diquark properties are 
consistent with those of other authors [241 22] • 

From Eq. m we can introduce an effective Lagrangian for the quark-diquark-nucleon 
coupling as in Ref. PH for the three-quark model 

PqDN = fs'^{ki)Osi’N(P)4>s(k2)A.s + fv^(ki)Oy(k2)—j^>N(P) ■ 4>Va(k2)A.y, (6) 

O s = 1, Oy = ( 7 " + z<9“/M) 75 , (7) 
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where ip, ipx, 4>s and cpy are the quark, nucleon and scalar and axialvector cliquark fields, 
respectively. The vertex functions Ag and Ay in Eq. © will be connected with the wave 
functions ip six, k±) and pv(x,k±) later. In this work we adopt the same form for A 5 and 
Ay. 

The Feynman rules of the diquark are presented in Table |T| 


Table 1: Feynman rules of the diquark. 
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The coupling of the virtual photon to a diquark is given as DUES! 

r 5 = -^Pi+PzT, (8) 

T Va0 = l{g<*p(Pl +P 2 Y - (! + K )9aijPl ~P2)p ~ (1 + k)$(jP2 “Pl)a} (9) 

with k the anomalous contribution of the magnetic moment of the axialvector cliquark, P 2 
the final momentum and pi the initial momentum of the cliquark. Here, we have used the’t 
Hooft-Feynman gauge parameter £ = 1, as in Keiner |25| , which ensures the zero charge of 
the neutron. 


3 Form Factors of the Nucleon 


The electromagnetic form factors of the nucleon are conventionally defined in terms of the 
electromagnetic current matrix elements, 


(N\'\r\NX) = -leuy(P') Y l Fi(q 2 ) + 7 - Q v [f, A^q 2 ) 

4777-tv 


MP) ( 10 ) 


with q = P' — P, P' and P the final and initial total momentum. 
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Figure 1: Feynman diagrams for the nucleon current in case a photon strikes the quark. The 
case of a photon striking the diquark is analogous. 


The form factors F)(g 2 ) can be determined from the J + matrix elements, 
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(ii) 


where Q\ = (puP^- + PLfPt)<i + ~ Zp^pfQL- If Q + = 0, the form factor is obtained in the 
so-called Drell-Yan-West frame. 


The Sachs form factors are 

G e ( 9 2 ) = U(9 2 ) + j^#-2(9 2 ), 

Gm(9 2 ) = Fi( 9 2 ) + C 2 ( 9 2 ). (12) 


The form factor problem is now reduced to calculating the matrix elements 

= f d*k A(fc')S+,A(fc 2 ) 

AA ' J (2n) 4 D(k' 2 )D(k)D(k 2 y y ’ 

as shown in Fig.|T| Here A(k 2 ) and A{k 2 ) are the vertex functions Ag or Ay given in Eq. 
and 


D{k) = k 2 — m\-\- is, 

D{k' 2 ) = k' 2 — m 2 -\-is — (P 1 — kj 2 — m 2 + is, 

D(k 2 ) = k 2 — m\ + is = (P — k) 2 — m 2 + is. (14) 
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The spinor part S xx , can be obtained from the Feynman rules in Tableland Fig. |T| If 
a photon strikes the quark, the quark is labeled as particle 2, with the spinor parts 

= -•(/|[*4]<VY'‘[fc]|*)> (15) 

= *S<jg^(f\6^\k' 2 )e,Ylh}0’v^\i) (16) 

in the numerator of Eq. ns for the current matrix element, where we use the notation 
[p] = 7 ■ p + m, as earlier for [P], and O = . 

If the virtual photon strikes the diquark, the diquark is labeled as particle 2 with the 
spinor parts 

sis = </IWIi)|rg. (17) 

siv = (/|0^(fc]O^|i>Tr[r j e' I r,]r^, (18) 

with rg, V^ a p from Eqs. (18191) . Compared with the three-quark model, only Fig. Ufa) and 
Q)d) of Ref. |E] remain in the diquark model because we regard the diquark as a pointlike 
particle so that such details of its internal structure are no longer relevant. 

From an explicit calculation we recognize the same isospin parts as in Kroll et. al. m 
For example, the isospin part of Eqs. (I15H18I) for the proton is e u , e u + 2e^, e u d, e u d + 2e uu 
respectively, where e u d, e^, e uu are the electric charges of the diquarks as given in Ref. 1221 

4 Zero Mode Analysis 

In order to search for zero modes, we now integrate the current matrix element over k~ in 
Eq. [mi in the valence and nonvalence regions corresponding to Fig. |T| (b) and (c) In frame 
where q + > 0 . 


7 


4.1 Integration in the Valence Region 


In the valence region, characterized by 0 < k + < P + , there is only one pole k~ = k~ n = 
(ml + k 2 L — i)/k + located in the lower half of the complex fc“-plane. Thus, we can integrate 
over k~ by Cauchy’s formula which yields 

+ _ [ d 4 k h(k' 2 )S+ y h(k 2 ) 

AV J (27t) 4 (k' 2 — ml + ie) (k' 2 — A' 2 + le)(k 2 — ml + ie) (k'l — A 2 + te) ( k 2 — m\ + te) 


= 2m / dx 


dk 2 . 


h(k! 2 )St>,,h{ki)\ k -. k - 


o 2(2v) l xxlx} (M] - - M\) 


(19) 


where Mg = Ml 2 = + *2^ Afi = + ^4, M(( 

U X2 X ’ U x X ' A X2 X ’ A 


2 i A' 2 


fc i +A‘ 


+ 




w ith x = /c + /P + , x 2 = 1 — x, x' = k + /P' + = x/a, x ' 2 = 1 — x', k' ± — k± — x'q^. 
Here we adopt Pj_ = 0, and Mf and M.\ are the mass of the initial and final states, i. e. the 
nucleon mass. We use a Pauli-Villars regulator A in A (k 2 ) = h(k 2 )/(kl — A 2 + ze), where A 
plays a role of momentum cutoff, to make the integration in Eq. (fTTH) and Eq. (J2H) below 
finite SEE. From Eq. m we see that, if all the poles arising from A(k 2 ), such as more 
Pauli-Villars fields, lie in the upper half of the complex V-plane, the explicit form of A (k 2 ) 
has no effect on the result in the valence region provided the definition as Eq. (|22j) are chosen 
to include all factor from A (k 2 ) to definition. A more explicit discussion of Pauli-Villars fields 
can be found in Ref. m in detail. 

We can connect the integration in the valence region obtained above to the form of 
Ref. jHH directly by using the spin sum 


2 m[p\ = E «a(p)«a(p). 


( 20 ) 


For example, 


= oc u / u(A:2)^=e ? 7 ,i ^Mu(A:2)uj 

y/ x 2 


( 21 ) 


ix' 2 x 2 yx 2 \/k 2 ' \jk!\ 

obviously is a current density between two vertices Ufu(k 2 ) and u(k 2 )ui i 


m 




















Eq. 21 Then we let 


1 t i \ 

V ^' 21 j 


( 22 ) 


where /i(A; 2 )/[:r 2 (m^ — M^)\ corresponds to A(fc 2 ) before integration. </?(x 2 , fc 2 j_) is just the 
usual radial wave function in eq. ©• So if we use a A(fc 2 ) without pole, our dehnition 
of <p(Y 2 , fc 2J _) is the same as the one used by Araujo et. al. [TH] , And only with such a 
dehnition do the x, X 2 x 2 that are left give an invariant phase-space volume element, clT = 
d 2 k±d 2 k2±S(k± + k2±)dxdx2/ (xx2)5(x + X 2 — 1) in the Drell-Yan-West frame fTBj . 

We find for the valence region results similar to Ref. [TH] , In the three-quark model the 
results with the integration as in this work are the same as the results in Ref. [HI because 
all terms with \p] are on-shell so that they can be evaluated as J2\U\(p)u\(p)- 

I 11 the q + = 0 frame, that is, the Drell-Yan-West frame, J xx , can be obtained by the limit 


= li m J™- 


a —>1 


AA'* 


(23) 


In the valence region corresponding to Fig. Q] (b), when a —> 1, then x’ — > x, x' 2 —■> x 2 , 
and so on. We then End that W is obtained from the wave function in Eq. (J2|) and the 
free quark current density ^j=e g 'y fJ, ^Y= or the diquark current in Eqs. dUEl) as usual. From 

V fc 2 

Eq. (fill we see that, if a pole of A (/c 2 ) arises from (k 2 — m 2 + xe), where m = m 2 , A, etc, it 
will lie in the upper half of the complex /c”-plane. Thus, the explicit form of A(/c 2 ) has no 
effect on the integration in the valence-quark region provided the diquark structure function 
h(k 2 ) is chosen appropriately, as will be discussed below. 


4.2 Integration in the Nonvalence Region 

The nonvalence region is characterized by P + < k + < P' + kinematically. We can recognize 
Fig. □ (c) as a pair creation process in the nonvalence region from the orientation of the 
momentum. Thus, the nonvalence contribution corresponds to a Z-diagram in the three- 
quark model if the virtual photon strikes a quark and corresponds to a double Z-diagram 
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if the virtual photon strikes the diquark. Thus dynamically, the nonvalence region can 
no longer be described by a wave function. In the Drcll-Yan-West frame there is no such 
contribution to the J + current component, so it must vanish when q + —> 0, that is, a —> 1. 
Otherwise there would be a so-called zero mode contribution. 

In the nonvalence region there are two poles located in the lower half of the complex 
A;”-plane when we integrate the current density over k~ by Cauchy’s formula to yield 


T+ — 
J XX' ~ 


d A k 


h(k' 2 )S+yh(k 2 ) 


(27t) 4 ( k' 2 — m 2 + ie)(k' 2 — A ' 2 + ie)(k 2 — m 2 + ie)(k 2 — A 2 + ie)(k 2 — m 2 + te) 

Kk' 2 )st x ,h{k 2 )\ k _ =p ._ k . 


= 2m / dx 


dki 


r{ 


f 2 on 


i 2(27r) 4 a;a;"V 2 A ' 2 - m\ (Mj - M' 0 2 )(q 2 - M m2m2 )(q 2 - M m2 \) 


h(k' 2 )Sl x ,h(k 2 )l 7 _ , 


Pf ~ k 2A 


( M} - M' 2 )(q 2 - M A , m2 )(q 2 - M A , K ) 
where x" = x 2 /{l— a), fc" = k±— x"q± } M m2m2 




(24) 


fc^+ra 2 , k^+m 2 itj k^+m^ , fej_ 2 +A 2 

-M-m 2A i ~ I “ ■ 


1 —a:'' 


1 —a;" 




fc " 2 +A '2 fe '; 2 + m 2 


1 — x n 


+ 


fc'; 2 +A ' 2 , fc'' 2 +A 2 


, ATa'a = —1- 


If there are more Pauli- Villars fields used, that is, A (k 2 ) = h^k^/Yli (*2 ~ A? + xe), we 
can use relation 


k 2 - A 2 + le , 


(25) 


( k 2 - A 2 + le) ( k 2 - A 2 + xe) A 2 - A 2 y A; 2 - A 2 + re 

to split the integration to the sum of integration as Eq. (12411 . 

Now we connect the integration in the nonvalence region contribution to the form of 
Ref. |E|. In the valence region we used the spin sum 2 m\p\ = J 2 x u x(p)ux(p)- However, in 
the nonvalence region we have Xi = x — 1 for initial part, so we should define k 2 as k — P , 
and not P — k as k 2 , that is, all k 2 should be replaced by —k 2 . For example, 


CM 

D qS 


= -*</1 


[k' 2 ] 


-k 


1X 2 X 2 

using the spin sum 


e q T 


V^2 


21 H) oc Ufu(k 2 )^Y= e q n f^~^= v (k 2 )ui (26) 


At A4 


2 m[-p\ = ~J2 v *(p) v x(p)- 
x 


(27) 
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Obviously, this is a current density ^-7=e q between two vertices Ufu(k 2 ) and v(k 2 )ui. 

V k 2 v fc, 2 

Hence the nonvalence region contribution is just the virtual pair creation process. 

Because h(k 2 ) functions in both regions are from the same one in Eq. (FTTH) we can still 
define h[k 2 ) as 

7 1 7—2- 2 - T72\ = v( x 2 ,k 2 ±)- (28) 

^/2(2vr) 3 ( m N - M$)x 2 {m 2 N - M\) 

However in Eq. (|2H) only (M‘j — M ' 0 2 ) and ( M 2 — M ’ 2 ,) occur for the final state. The vertex 
in the initial state is the so-called no-wave function vertex which, via Eq. 0 , has a spin part 
different from the nucleon wave function but still a radial nucleon wave function. This vertex 
was also studied in Refs. ezi 12a mi- In the work by Choi et. ah they use BS equation to 
connect the nonwave function vertex with the normal vertex with a kernel. In our work, we 
do not make explicit calculation of nonvalence contribution, so we use the results without 
using BS equation directly. 

4.3 Recognition of Zero Modes 

Now we proceed to look for a zero mode of the nucleon, that is, the spin-1/2 case. The zero 
mode contribution is defined as in Ref. EQ 

= li 11 ) J\\'nv (29) 

No singularities, hence zero modes, arise at x — 0. To isolate the zero mode contribution at 

x = 1 or x = a, we first make a transformation of the longitudinal integration variable, 

t+ r r*j x(a - x)[x" {I - x")] 2 / , Mc+ un u ! 

Jix'z.M. ~ h m / dx -- h(k 2 )Sl x ,h{k 2 ... 

«-*lj 1 XX Z X 2 

= lim f dx ^— ^-h(k' 2 )SZ x ,h(k2)[-] 

a—>l J 1 [1 — ay 

= lim [ dza(a — 1)(1 — z) 2 h(k'yS^ x ,h(k 2 )[...} (30) 

a-*l Jo 

where x = 1 — z( 1 — a) and x" = (1 — x)/{l — a) = z. So, both a — x and 1 — x can give 
a factor 1 — a. Hereafter we identify a — x and 1 — x, when we consider the zero mode 
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contribution. Our analysis at the endpoint 1 — x will be based on counting powers of 1 — x 
in the integrand of Eq. ©. 

Now we take 

1 A 2 

^ X2,k2± a /2(2tt ) 3 (m 2 N - M$)x 2 (m 2 N - M \)' 31 

that is, h(k 2 ) = A 2 and h(k 2 ) = A /2 , which is just the form used by Choi et. al |21j as a 
starting point of their analysis, because this is the simplest choice and it means that there 
is no complicated internal structure in the vertex. Then the key to recognizing a zero-mode 
contribution is the spinor part S xx , in Eq. (1101) . Using Table IV in Ref. ra, it is easy to 
count (1 — x) n factors in S xx ,, with (1 — x) coming from x' 2 or x 2 . Every [k 2on ] or [k 2on ] in 
S xx , leads to a term with the factor (1 — x) _1 , where k 2on or k 2on are the momenta before 
and after the photon strikes the quark or diquark. Adjacent [, k 2on ] and [ k 2on ] terms (that is, 
without another [p on \ between them) give a single factor (1 — x) because Af k or L in Eq. 
(All) in Ref. [T2J will give a (1 — x) factor when q —> 0. The y + and y _ between [k 2on \ and 
[k 2on \ give a factor (1 — x) and (1 — a:) -1 , respectively, and y + between [ k 2on ] and [ k 2on ] gives 
a factor (1 — x). In contrast, y + and y~ in other cases and other y’s, such as y5, y L , do not 
give rise to terms with factors (1 — x) n . 

Applying this power counting technique to the following example, when the spinor part 
is given by m 

s h = ?M[/4]y + (l-y 5 )[fc 2 p]y 5 e* ’ r ] 

= Tr [[ k 2onh + ( l ~ lb)[ k 2on]{konY ' Ty 5 ] 

+ Tr [W2on)l + ( l ~ lb)[k2on]l + e* ■ Ty 5 }^(k~ ~ kon), (32) 

where T M = y M — {k + k')^ l /D and e Q is a polarization vector, then in both lines after the 
second equality sign, the (1 — x)~ 2 from [k 2on ] and [k 2on ] is canceled by y + (l — ys) and the 
adjacent [k 2on \ and [ k 2on ] terms. For the e* ■ y part of e* ■ T in the first line, there will be 
(1 + x) and (1 — x)° from e* - y + y5 (or e* - y“y5) and e* • yy 5 respectively. However e* ■ {k + k') 
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in the D term will give a (1 — x ) _1 term. For the e* ■ 7 part of e* ■ T in the second line, the 
7 + from e* ■ V disappears because of the other 7 + using 7 +2 = 7 + . The trace will let [k' 2on ] 
and [k 2 on] become adjacent again. So, the (1 — x)^ 1 from \{k~ — k on ) will be canceled. Then 
there will be (1 — x) and (1 — x)° from 7 + e* • 775 and 7 + e* + 7 _ 7 s, respectively. Thus, a zero 
mode has to come from D in in (1 — x)° order. 

We now use the same method to analyze S q s, S q v, Sds, Sdv of Eqs. (11511181) . 

For S q s the analog of Eq. (H?T) is [k' 2 \e q 7 ^[A; 2 ] from Eq. (TTlol) in which the (1 — x )~ 2 from 
W 2 on] and [k 2on \ is canceled by the y + matrix element and the adjacent [k' 2on ] and [ k 2o n ] terms. 
The same applies to Sqy , where the analog of Eq. (B21) is Oy[k' 2 \eq'-f ll [k 2 }Oy from Eq. (ITU . 
For Sds hr Eq. m the instantaneous part of the term [k], which comes from the quark 
propagator, is q + (M 2 — Ad' 2 )/2P' + , which will give a (1 — x)~ l factor. This (1 — x ) _1 factor 
will be canceled by that from ( k 2 + k' 2 ) + . Since there is already a factor 1 — a in Eq. 1301 
there will be no zero mode from S q s, S q y and Sds- 

Now we consider the spinor part Sdv whose analog of Eq. (ET2|) is Oy^[k]Oy^Tr[Tie q Tj] 
of Eq. (TTCrl) with from Eq. (J31). As far as g a / 3 (pi +P 2 ) fl hr Y ,J a g is concerned, ( p\ +p 2 Y can 
cancel the (1 — x)~ l factor from [k] as in Sds- However, for (1 + — gpq a ), the factor 

(1 — x)~ l can not be canceled. Hence, if we remember that h{k '. 2 ) and h{k 2 ) are constant, as 
in Choi et. al. [2Q, there is a zero mode contribution in (1 — x)° order from Sdv- 

From the above analysis, with a constant h(k 2 ), there will be zero mode. Obviously, if we 
choose ip(x 2 , k 2 ±) to let h{k 2 ) have a pole, that is, have a factor (1 — x) n with n < 0 , there 
will be a higher order zero mode while if h{k 2 ) has a factor (1 — x) n with n > 0, which could 
reduce h{k 2 ) after the transformation in Eq. (1HU1) to zero, the zero modes will disappear. 

In the current matrix element, Eq. (El, the term [k] comes from the quark propagator 
whose instantaneous part is 7 + (M 2 — Ad' 2 )/2P' + . Hence the corresponding instantaneous 
spinor part 

Stv,n„ = (/l<5^W.„ ai Ot^|i)Tr[r ie ,r,]r7 
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T*t T j 


= i{N\-^-^\N)Tr[T i e q T j ]{M 2 - M ,2 )/2P' + 
• ma'7 5 (7“ + P'7^)7 + (7 /3 + PP/M^ux 


■ {gapijPl + P 2 ) + + (1 + «)#+ ( k 2 ~ fe)/9 ~ (1 + {K ~ h)a} 

3 

~ *(1 + ^(iVl^^l^rr^^M 2 - M' 2 )/2M 
• -u A /[7 + 7-?-7-97 + -7 + ? 2 / m ] m a 


= C(l +k)(M 2 - M' 2 ) 


P+q 2 /M 2 , A'= 1/2, A = 1/2 
2P + q L /M, A' = 1/2, A = -1/2 


(33) 


where C = ?(A r |^-^|A^)Tr[r i e 9 r J ], M' 2 = M/ 2 or M' 2 ,. After ~, we have omitted all higher 
order terms of (1 — a) because they will vanish when q + —> 0. Hence there is clearly a zero 
mode in Sdv , when g + —> 0. 


4.4 Numerical Estimate of the Zero Mode 

To give an estimate of the zero mode, in Fig |21 we present the valence and nonvalence 
contribution of Jdv(Q 2 )/^dv( 0 ) i 11 the Q + —► 0 limit with h(k' 2 ) and h(fc 2 ) taken as constants 
using Eqs. CES El ESD- The valence contribution can be calculated using Eqs. (fiBi m 
directly. The nonvalence contribution can be obtained from Eqs. (in ei and the terms 
involving a zero mode have been given in Eq. Here we use a mass rn q = 0.43 GeV 

for the u (d) quark as in Refs. [21, 2T] and a diquark mass of mr, = 0.6 GeV. We also 
adopt A = A' = 1.5 GeV. For both cases, the nonvalence contribution is visible. In the 
case A = 1/2, A' = 1/2 in Fig. 2(a), the nonvalence contribution increases with Q 2 and will 
cross the valence contribution at about — q 2 = 4 GeV 2 . The non valence contributions will 
decrease slowly after about Q 2 = 8 GeV 2 . Hence the zero mode is far from negligible. 

This power counting method we used is not required in the spin-1 meson case because 
the trace can be evaluated. However in the spin-1/2 case, this is difficult to do especially in 
the three quark model. 
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Figure 2: JpviQ 2 ) th e valence and nonvalence regions divided by <//>y(0) l 11 th e valence 
region, (a) is for the case A = 1/2,A 7 = 1/2 and (b) for the case A = 1/2,A 7 = —1/2. The 
full line is for the valence contribution, the dashed line for the nonvalence contribution. 

5 Numerical Results 

In Sec. 4, we have found that if the wave functions are chosen to let h{k 2 ) and h(k 2 ) be 
constants, there will be a zero mode arising when q + —> 0, i. e. the Drcll-Yan-West frame 
is used, which indicates that the nonvalence contribution must be considered. In the meson 
case, there are many works calculating a nonvalence contribution [23 ESI EH]. It would also 
be interesting to study the nonvalence contribution in the timelike region where g 2 > 0 so 
that we must use a frame with q + > 0. However in this paper, we focus on the spacelike 
region where q 2 < 0 so that Drell-Yan-West frame can be used. I 11 case h{k 2 ) and h(k' 2 ) are 
not constant, from the analysis in Sec. 4 we can choose these wave functions to remove zero 
modes. 

We calculate numerical results in the Drcll-Yan-West frame defined by q + = 0. In this 
work we use for the wave functions <ps(x, k±) and <pv(x,k±) in Eq. 4 both the Gaussian 
parametrization 

<p{x,k±) — N exp(—M q/8/3 2 ) , (34) 
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and the negative power parametrization indicated by the form factor of Ref fSIJ • 

ip(x, k±) — N {to 2 + Mq ) -3,5 , (35) 

where N is a normalization constant. Inserting the Eq. (ETH) to Eqs. (EU EH) from Eq. (EDD, 
we can find that the Gaussian form removes any zero mode because of Mq ~ {x — 1 ) _1 (Here 
we note that in the nonvalence region x^i = x — 1 > 0 and X 2 f = oi — x > 0 as given in 
fig- Kb), which ensures Mq > 0 when q + —> 0). The negative power shape wave functions 
will also remove the zero mode because h(k 2 ) and h(k' 2 ) determined by this wave function 
will generate a (x — l) 5 factor in Eq. (|3UJ) . 

The mass of the quark and diquark are chosen as m ^/3 and 2mjv/3, respectively, and 
k = 1.6 as in Ref. j25j (in fact, the results are not sensitive to k). From our calculations we 
End the results are more sensitive to the coupling constants fs and fy. In the works using 
direct Melosh rotation |111131 j . the proton charge can be obtained automatically. However, as 
discussed in Sec 2 , we let the scalar diquark coupling constant fs determined by G P C (<S) = 1. 
In Refs. [HH I32j , they also normalized the wave function to charge one of proton. The 
free parameters left are fy , f3 for the Gaussian type and to for negative power type wave 
functions. The zero charge of the neutron is guaranteed by the formula we used, that is, it 
is not fitted. 

The parameters characterizing various fits are given in Table II. 

Table 2: Parameters for Gaussian type wave function (Gau) and for negative power one 

(NP). 



m q (MeV ) 

rriD(MeV) 

K 

/ 3{MeV) 

to(MeV) 

Gau 

313 

626 

1.6 

225 

— 

NP 

313 

626 

1.6 

— 

100 


We present the charge and magnetic radii and magnetic moment of the nucleon with 
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different fy in Table III. 


Table 3: Charge and magnetic radii and magnetic moments for different values fy. The data 
are from Refs [331 E3 j SI EH)- 



fv(MeV) 

r c(f m ) 

r P M U m ) 

rp n (/m 2 ) 

r n M (fm) 


Mn(/^Ar) 

Gau I 

80 

0.804 

0.834 

-0.262 

0.861 

2.225 

-1.009 

NP I 

80 

0.821 

0.855 

-0.274 

0.882 

2.228 

-1.013 

Gau II 

100 

0.821 

0.820 

-0.284 

0.850 

2.427 

-1.075 

NP II 

100 

0.836 

0.837 

-0.295 

0.868 

2.428 

-1.078 

Gau III 

130 

0.854 

0.796 

-0.325 

0.832 

2.814 

-1.201 

NP III 

130 

0.872 

0.816 

-0.339 

0.852 

2.811 

-1.203 

Exp. 

— 

0.870(8) 

0.855(35) 

-0.116(2) 

0.873(11) 

2.793 

-1.913 


As in the work of Kroll et ah mt the results for the proton are better than those for the 
neutron. The charge and magnetic radii and magnetic moment are close to the data with 
fy = 130 MeV. The magnetic moment of the neutron is too small, while the charge radius 
is too large, as is the case with other authors. 

The electromagnetic form factors of the nucleon are shown in Fig. [31 Since the experi¬ 
mental magnetic form factors and the electric form factor of the proton are well described 
by a dipole fit 

G r E (Q 2 ) = G V M {Q 2 ) / jj, p = G" M (Q 2 )/ ft , = F d (Q 2 ) = (l + IT) (36) 

V m Dj 

with ml = 0.71 GeV 2 . 

In Fig. [31 we present G P E {Q 2 ), G P M (Q 2 ) / /i p and G 1 f I (Q 2 ) / fj, n divided by F E (Q 2 ) and the 
charge form factor of the neutron, G E (Q 2 ) //i p , with fy = 100 MeV. The magnetic form 
factors of the proton and neutron are close to the dipole form with both wave functions 
while G p m (Q 2 )/ /j, n is lower than the dipole form, which is consistent with Jefferson Lab 


17 













Figure 3: Form factors of the nucleon, (a) and (b) are for the Gaussian type and negative 
power type wave functions with fy = 100 MeV, respectively. 

data maim eh. The magnetic form factor of the neutron is a bit too high. 

In Fig. flU), we present the ratio iiG p e /G p m with different fy compared with recent data 
from Jefferson Lab. With all fy , the ratio will decrease in the higher Q 2 region and the falloff 
is too fast. There are similar results for the three-quark model in Ref J21J. The difference 
between the data and experiments may arise from the simplistic wave functions we use and 
omission of the intrinsic form factors of quarks and diquarks. In Fig. Q we also find that the 
momentum transfer Q 2 , where the electric form factor of the proton crosses zero, increases 
with increasing fy. 


6 Summary and Discussion 


We have given a relativistic diquark model of the nucleon using the KW method. We calcu¬ 
late the electromagnetic form factors of the nucleon taken as a scalar and axialvector diquark 
bound state. In this work we do not consider the effect of scalar-vector transitions and pion or 
gluon exchange effects, which may improve the numerical results. The charge and magnetic 
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1.2 



Figure 4: The ratio /iG v E jG v M . The thick and thin lines are for the Gaussian type and 
negative power type wave functions with f v = 80 MeV (full), 100 MeV (dashed), 130 MeV 
(dotted) respectively. Experimental data are from Refs. [STJ 3B, E9|. 

radii and magnetic moment of the proton can be reproduced, while the magnetic moment 
of the neutron is too small. The shape of the form factors G P M (Q 2 ) //i p and G^(Q 2 )//i n can 
be reproduced also. The zero charge form factor of neutron at Q 2 = 0 (neutron charge) is 
guaranteed by the formalism we used. 

There will be a zero mode in electromagnetic form factors of the nucleon arising from the 
instantaneous part of the quark propagator if we adopt ip(x 2l k 2E ) to let h(k 2 ) have (1 — x) n 
with n < 0. In our numerical calculation, we find that the zero mode cannot be neglected. 
Hence the quark-diquark-nucleon vertex functions must be chosen to remove the zero mode 
to ensure the correct results in the Drell-Yan-West frame. 
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